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Another bar in the Bulge. 
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Abstract. A map of the projected density of the old stellar population of the Galactic Bulge region is reconstructed 
using 2MASS data. By making a combination of the H and K photometric bands , it is possible to overcome the 
effect of reddening, and thus penetrate the inner structure of the Galactic Bulge. The main structure in the map 
corresponds to the well documented peanut shaped bar which is formed by the inner parts of the Galactic disk as 
a result of dynamical instabilities. As suggested by numerical simulations, the projected Z profile of the bar, has 
an almost exponential shape. After subtracting the exponential profile associated with the bar, a large residual 
appear near the Galactic Center. This residual is elongated and asymmetrical, which suggest a bar structure. 
Thus we arrive at the conclusion that in addition to the main bar a smaller bar with a different orientation may 
exist in the central region of the Milky Way. This finding makes the Milky Way very similar to a large number of 
barred spiral Galaxies which show as well a smaller bar in their central regions. 
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The recent release by the 2MASS collaboration of a cata- 
logue of point sources covering a large fraction of the sky 
is of great interest for the study of the Galactic Structure. 
Since the 2MASS infrared magnitudes by 2MASS arc in- 
sensitive to extinction, these data are an ideal tool to 
probe the structure of the Galaxy within a few degrees 
from its center. The study of the Galactic Bulge in low 
extinction windows situated at larger distances from the 
center revealed that the Bulge was flattened and asymmet- 
ric. This elongated bulge with marked tri— axial structure 
is usually called the "bar" . The major axis of the bar 
is pointing towards the first Galactic quadrant. The pres- 
ence of the bar is very noticeable in the stellar component. 
Stanek (1994) detected an offset in magnitude for clump 
giants on different sides of the bar. A Similar asymmetry 
is visible in the COBE image of the Bulge region taken in 
the near infrared (Dwek et al, 1995). But the bar is also 
quite visible in the pattern speed of the gas (Binney et al, 
1991). The nature of the Bar/Bulge system and its for- 
mation process are not completely clear. Spiral Galaxies 
like the Milky way have exponential bulges (Carollo et 
al. 2001, Mollenhoff, C, Heidt, J., 2000, Van der Marel, 
Roeland P., 2000), since exponential bulges can be formed 
easily has a result of disk instability in the central regions, 
we may infer that the exponential bulges are formed by a 
dynamical perturbation of the central disk. This bar struc- 
ture seems to dominate the central region of our Galaxy, 
but a fundamental question is to know whether other corn- 
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ponents do exist at smaller scales. One may expect to find 
some parts of the bulge that were not formed dynami- 
cally from disk stars, but from a distinct stellar popula- 
tion. This population should be dominant in the central 
degree, and should be detectable using the 2MASS data. 

2. Data analysis 

2.1. Basic processing 

The 2MASS catalogues corresponding to the range \bu\ < 
10 and |/jrj| < 15 were obtained from the 2MASS public re- 
lease. In this release the coverage of this coordinate range 
is not complete (~ 60%), but is sufficient for studies of 
Galactic structure. The first step is to select the sources 
with sufficient quality (read_flag>l). Using this criterion 
there are about 30 millions stars in the region of interest. 
Within this area there are a number of regions which are 
affected by the presence of bright stars and their diffrac- 
tion spikes. It is easy to identify these areas since the star 
counts in the neighborhood are much lower than average. 
By making star counts in a box of 0.15 sq. degree all over 
the frame, we obtain an image where the regions occupied 
by bright stars appear as dark patches. In case the counts 
in a box are less than 10, the relevant pixel in the im- 
age is considered to belong to a dark patch. Even in the 
lower density regions, the mean star counts are about 10 
times larger, thus our cut-off does not induce any artifact. 
To remove any contamination from the dark spots to the 
nearby pixels, all the pixels belonging to a 3x3 mesh were 
also flagged. 
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2.2. Canceling the effect of extinction. 

By combining 2 photometric bands it is easy to derive a 
magnitude that is independent of reddening. For instance 
using H and K, we can derive the apparent magnitude m e : 

m e = K- AK (H-K) 
Ah - Ak 

Following Rieke and Lebofsky (1985) we take: 

4^ = 0.112 and 4^ = 0.175 

A v A v 

Leading finally to: 

m e = K-C K {H-K) with : C K = 1.77 (1) 

Note that since the reddening law in the infrared is 
unique (Cardelli et al. 1989), the definition given in eq (1) 
is valid for any line of sight, and not only for the par- 
ticular area probed by Rieke & Lebofski (1985). Thus in 
this band the cumulative counts of the stars brighter than 
some cut— off in magnitude Mc will be independent of the 
reddening (provided that Mc is brighter than the obser- 
vational limiting magnitude). Considering that the limit- 
ing magnitude in the H band is about 16 and that the 
maximum extinction is about Ay ~ 30, we infer that the 
cut— off should be brighter than H = 11. Since the color 
of the tip of the upper giant branch is H — K ~ 0.5 we 
obtain finally in the m e band, Mc — 9.6. For security we 
will adopt the very conservative bright upper cut— off in 
magnitude Mc — 9. 

2.3. Reconstructing density maps of the bulge region. 

There is one major difficulty in producing a map of the 
Bulge region: the coverage is not complete and addition- 
ally there are holes due to the bright stars. However pro- 
vided that we assume that the Galaxy is symmetrical 
about its plane, it is possible to fill most of the gaps. In 
order to smooth and fill the smaller data gaps which re- 
main, we use the fact that the density profile at constant 
longitude is almost exponential (see Fig. 1). This property 
suggest that the profile can be represented by a polyno- 
mial function. Numerical experiments shows that it is not 
significant to increase the degree of the polynomial be- 
yond 5 to represent the data. To increase the numerical 
stability of the fit we will use a strip of 9 columns centered 
around the column of interest. This procedure is carried 
out for each column in the image. A filtered image is con- 
structed by replacing each column in the original image by 
the polynomial solution. Some example of polynomial fit- 
ting of the columns are given in Fig. 1. This procedure has 
a good ability to fill or extrapolate small data gaps, and 
has an excellent numerical stability. Once we have recon- 
structed this image of the star counts in the Bulge region, 
we apply a final wavelet smoothing procedure in order to 
balance the smoothing at all scales. This final smoothing 
is interesting because the polynomial reconstruction is one 




2 4 6 8 

b 




2 4 6 8 

b 



Fig. 1. Here we present the polynomial smoothing of 2 
sections of the projected density taken at constant longi- 
tude. Note that the profile deviates from an exponential 
near the Galactic Center. 

dimensional, and thus is biased in one direction, on a par- 
ticular scale. The Wavelet decomposition is obtained by 
applying iteratively the Spline filter to the image and the 
smoothed images (Starck & Murtagh, 1994). To estimate 
the statistical cuts to apply in the wavelet decomposition 
we generate Monte- Carlo images with counts approaching 
our own image. In the final reconstruction of the smoothed 
image we use 4 a cut— off. One final concern is the pos- 
sible effect of small uncertainties in the determination of 
Ck (typically a few %). This can be investigated by re- 
constructing the density with small variations of Ck- The 
comparison with the initial map shows that the density 
variations induced by such changes in the value of Ck are 
about the amplitude of the noise, and thus do not affect 
the final result. 

2.4. Analysis of the projected density. 

We have already noticed that the Z profile of the projected 
density is almost exponential. This exponential profile is 
also present in numerical simulations of peanut shaped 
bars. Combes et al. (1990) showed that a disk with a small 
bulge near its center forms a peanut shaped bar with a 
nearly exponential projected Z profile (see in particular 
Fig. 4 in Combes 1990). Thus if we subtract the expo- 
nential contribution which corresponds to the bar— disk 
system, the remaining density may reveal another com- 
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Fig. 2. Map of the Galactic Bar region reconstructed using a polynomial fitting method and wavelet smoothing. Note 
that the star counts are systematically higher at positive longitude for \b\ > 2. Contours values: (max,min)=(60000,400) 
stars/sq. deg. 
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Fig. 3. Contours of the residual density after subtracting the main component. Note that the asymmetry is in the 
opposite direction of the large scale Galactic Bar. Contours values: (max,min)=(38000,400) stars/sq. deg. 



ponent. The contribution of the bar will be estimated by 
fitting an exponential to each column of the image (which 
corresponds to the projected Z profile). This procedure is 
more flexible than trying to subtract a given bar model. 
There are still many uncertainties concerning the struc- 
ture of the Galactic bar, thus the subtraction of a bar 
model may give ambiguous results. By fitting an exponen- 
tial profile, we make no particular assumption about the 
shape of the bar, other than a general assumption on the Z 



profile at equilibrium which is justified by numerical sim- 
ulations. To implement the fit of the exponential profile, 
we perform a robust fitting of a straight line to the log of 
density (by minimizing the sum of absolute deviations). 
Once the exponential contribution has been subtracted, 
a very significant residual appears in the central region 
(R < 2°). The contours of this residual are smooth and 
elongated along the Galactic plane. This component shows 
also a very significant asymmetry in longitude. The am- 
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plitudc of this asymmetry is close to 15 %, which is about 
7 sigmas according to Poisson statistics. There are also 
some residuals along the Galactic plane in general. But 
their amplitude is about 10 times smaller, they do not 
have smooth structure, and their scale is much smaller. 
These residuals are probably due to the presence of young 
stars in the HII regions. 

3. Discussion 

There are some intrinsic difficulties in the interpretation of 
the former results, we see a residual component after sub- 
tracting a main component. But this result is based on the 
assumption that the bar has a projected density profile in 
the Z direction which is exponential. Even if this assump- 
tion is supported by numerical simulations, it might just 
be that the projected density of the bar is not exponen- 
tial. But in this case, why do we observe an asymmetry in 
an opposite direction to the bar ? Two possible sources of 
bias are: the residual extinction, and the steepness of the 
density profile near the center (Blitz & Spergel 1991). In 
principle the extinction should have no effect in the the 
band, but it is possible that for observational reasons the 
limiting magnitude in 2MASS is somewhat brighter than 
expected, which would result in an indirect extinction ef- 
fect. These biases can be investigated by using numerical 
simulations. To build a numerical model we need to inte- 
grate the convolution product of the luminosity function 
<t> with the density distribution p. The integration domain 
in the space of the magnitudes will be modulated by the 
extinction Ay. For the luminosity function we will adopt 
the model of Wainscoat et al (1992). The density distri- 
bution will be built using a truncated exponential disk 
(Lopez— Corredoira et al. 2001) and a triaxial bar model 
with a power law profile. This bar model has an incli- 
nation of 20 deg. with respect to the line of sight, and 
axis ratio: ^ ~ 3 and ^ ~ 0.7 (Dwek et al. 1995). And 
finally the extinction map of the whole area was built 
from our 2MASS data by using a method presented by 
Schulthcis et al. (2000). Let's start with the case of resid- 
ual extinction effects: the counts were generated using the 
aforementioned procedure, Poisson noise was simulated, 
and finally the whole process of polynomial reconstruc- 
tion, smoothing and exponential subtraction was applied. 
This procedure was undertaken for different observational 
limiting magnitudes, starting from our default value (no 
extinction cut— off). To summarize the results, the longi- 
tude profile of the residual density has been represented 
for the different limiting magnitude (Fig. 4). The compar- 
ison of these results with the observational profile shows 
unambiguously that the effect of extinction is in the oppo- 
site direction to the observed asymmetry. Furthermore, it 
is interesting to notice the steep edge of the observational 
profile, which is similar on both side of the diagram. The 
profile variations produced by extinction are not similar on 
both edges, and are hardly so steep. Concerning the effect 
of the steepness of the inner bar profile, even if a steep bar 
profile, may produce an asymmetry opposite to the bar, 
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Fig. 4. The marginal distribution in longitude of simu- 
lated profiles for different limiting magnitude (default: 
thin line, 1 mag. brighter than default: dashed line, 1.5 
mag. brighter: dotted dashed line). The last profile is the 
observational profile (thick line). All profiles have been 
normalized so that the sum of the profile is unity. 

it is not possible reproduce the amplitude of the asymme- 
try even if the density profile behaves like: R~ 7 .A steeper 
density profile appears unphysical. These results are not 
very dependent upon the luminosity function. Power law 
luminosity functions gives the same result. The best expla- 
nation is to consider the presence of another component 
in the inner Galaxy. A small bar with a steeply dropping 
density near its edge can reproduce the observed asymme- 
try. This finding is not surprising, the discovery of small 
substructure in the central region of barred spiral Galaxies 
is very common (Erwin & Sparke 1999, Fiedli et al. 1996, 
Friedli 1996, Shaw et al. 1995). This small bar structure 
was not found in previous studies for 2 fundamental rea- 
sons, the former data set lacked either the depth or the 
resolution. The 2MASS H band is deep and insensitive to 
extinction, and furthermore the 30 million stars available 
result in an excellent spatial resolution. 
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